SPhT-T04/034 



Color Glass Condensate and BFKL dynamics 
in deep inelastic scattering at small x 

Kazunori Itakura 

Service de Physique Thcorique, CEA Saclay, F-91191 Gif-sur-Yvette, France 

o 
o . 

£SJ , Abstract. The proton structure function F 2 (:e,<2 2 ) for x < 10 2 and 0.045 < Q 2 < 

$_j ' 45GeV 2 , measured in the deep inelastic scattering at HERA, can be well described 

within the framework of the Color Glass Condensate. 



1. Introduction 

(N . 

There has been a surge of theoretical and experimental interest in the "Color Glass 
' Condensate (CGC)" which appears in the new perturbative regime of QCD relevant 
for high energy scattering pQ. This new state is characterized by high density gluons 
whose longitudinal momenta are very small compared to the total momentum of the 
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parent hadron (such gluons are called "small-x" gluons since the ratio of the momenta 
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is denoted as x). The gluon density is typically as large as and cannot be too 

large ( "saturated" ) so that the unitarity of physical cross sections is ensured. 

Recent rapid theoretical progress in understanding the physics of CGC is triggered 
by the experiments currently investigated at HERA (DESY) and RHIC (BNL). The 
relevant process at HERA is the deep inelastic scattering (DIS) of an electron off a 
proton, while at RHIC it is more complicated Au-Au or deuteron-Au collision. These two 
apparently different experiments (different in complexity and energy) are nevertheless 
closely related to each other in the context of CGC through the "universality" of the 
hadron/nucleus wavefunctions or the gluon distributions. Here we mean differently 
by the word "universality" than in the usual sense. Namely, in the saturated regime 
at high energy, the gluon distributions of the proton and nucleus are both described 
by the same function of the ratio of the transverse momentum of gluons k± to the 
saturation momentum Q s , which is the (inverse of) typical transverse size of gluons 
(for more explanations, see for example, Refs. [21 |Hj). Therefore, one will be able 
to 'translate' the HERA physics for protons into the RHIC physics for gold nuclei. 
Furthermore, the saturation scales in two experiments are accidentally of the same 
order because of its particular dependencies upon x and the atomic number A, i.e., 
Q 2 s (A,x) oc A 1 / 3 (l/x) A ~ (A/x) - 3 . Indeed, at HERA, A = 1 and x ~ 10~ 4 while at 
RHIC, A ~ 200 and typically x ~ 10~ 2 . This fact also strengthens the importance of 
understanding the HERA physics in relation to the RHIC physics. In this talk, I will 
show that the recent HERA data [3] at small-x with not so large Q 2 is consistent with 
the current picture of the CGC (for more detail, see Ref. which, according to the 
discussion above, suggests the importance of CGC at RHIC. 
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2. DIS at small x and previous fits with saturation models 

DIS at small-x looks simple in the "dipole picture" which leads to an intuitively 
understandable factorization formula. The scattering between the virtual photon 7* 
(emitted from the projectile electron) and the proton is seen as the dissociation of 7* 
into a quark-antiquark pair (the "color dipole" ) followed by the interaction of this dipole 
with the color fields in the proton. Then one can write the F 2 structure function as 

F 2 (x,Q 2 ) = I dzd 2 r\y TiL (z,r,Q 2 )\ 2 a dipole (x,r), (1) 

47T Qt em T L 

where, ^t,l are the light-cone wavefunctions of 7* with transverse, or longitudinal, 
polarization, and adi po ie(x,r) is the cross-section for dipole-proton scattering (for a 
dipole of transverse size r), containing all the information about hadronic interactions 
such as the unitarization or saturation effects. 

A simple parametrization for Odipoic^, r ) which has qualitatively plausible behaviors 
like color transparency and saturation effects was first proposed by Golec-Biernat and 
Wusthoff (GBW) jHj. They used a very simple function Odipoie^, r ) = °"o (1 — e -1 " 2 ' 9 ^ 1 ^ 4 ) 
with only three parameters <7o (a hadronic cross-section), Xq and A for the saturation 
momentum Q 2 (x) = (xq/x) x GeV 2 , and managed to provide rather good fits to the 
(old) HERA data for x < 10~ 2 and all Q 2 . This success was quite impressive by itself 
since it suggested the relevance of saturation physics in the HERA data, but at the 
same time required more serious theoretical effort towards understanding the HERA 
data with the saturation picture better rooted in QCD. In fact, there is no kinematical 
regime in which the GBW model can be (strictly) justified from QCD, and the GBW 
model must be improved with the information of QCD, or replaced by other QCD-based 
parametrization. So far, there are several attempts to improve the GBW model (ZHSj, 
but they mostly focused upon improving the behavior of the fit at large Q 2 , by including 
DGLAP dynamics. On the other hand, we know that the BFKL dynamics, rather than 
DGLAP, should be the right physics in the transition regime towards saturation. This 
BFKL physics was not addressed so far, and we will focus on the regime with not too 
large Q 2 where the BFKL and saturation physics should be more relevant, and will 
present a new analysis of the HERA data, which is rather orthogonal to the previous 
attempts to improve the GBW model. 

3. The CGC fit 

We restrict ourselves to the kinematical range where one expects important high density 
effects — namely, x < 10~ 2 and Q 2 < 50 GeV 2 — , and show that the data in this range 
are consistent with our present understanding of CGC (BFKL evolution and saturation). 
The upper limit on Q 2 has been chosen large enough to include a significant number of 
"perturbative" data points, but low enough to justify the emphasis on BFKL, rather 
than DGLAP, evolution. Within this kinematical range, we shall provide a reasonable 
fit (which we call the "CGC fit") to the new HERA data for F 2 based on a simple, 
analytic, formula for the dipole scattering amplitude. 
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Figure 1. The F 2 structure function as a function of x in bins of Q 2 for Q 2 < 15 GeV 2 
(left) and for Q 2 > 15 GeV 2 (right). The experimental points are the latest data from 
the HI and ZEUS collaborations gj. The full line shows the result of the CGC fit 
with 7V = 0.7 to the ZEUS data for x < 1CT 2 and Q 2 < 45 GeV 2 . The dashed line 
shows the predictions of the pure BFKL part of the fit (no saturation). In the bins 
with Q 2 > 60 GeV 2 , the CGC fit is extrapolated outside of the range of the fit. 



The dipole cross-section in the CGC fit reads <7 dipole (a;, r) = 2nR 2 M(rQ s1 Y) 1 
where the dipole scattering amplitude M(rQ s , Y) is constructed by smooth interpolation 
between two limiting solutions to the non-linear evolution equations in QCD [01 ED]: the 
solution to the BFKL equation with saturation boundary jTTJ [T2J H3] for small dipole 
sizes, r <C 1/Q s (x), and the Levin- Tuchin law [El Ej for larger dipoles, r 3> 1/Q s (x). 

AT(rQ s , Y) = A/- (^j for rQ s < 2, 

M(rQ Sl Y) = 1 — e - aln2 ( 6K W for rQ s > 2, (2) 

where Y = ln(l/x), Q s = Q s (x) = (xq/x) x ^ 2 GeV, and we have defined Q s in such 
a way that M{rQ Sl Y) = Ao for rQ s = 2. The coefficients a and b are determined 
uniquely from the continuity of J\f(rQ s , Y) at rQ s = 2. In the first line, 7 S = 0.63 (or 
more strictly, 1 — j s ) is the anomalous dimension, and k = x" (is) / x' (is) — 9.9 is the 
diffusion coefficient. The anomalous dimension gives the geometric scaling fH>\ ITT| IT2"]. 
while the second "diffusion" term in the power (the term depending upon k) brings in 
scaling violations. The overall factor Ao is ambiguous, reflecting an ambiguity in the 
definition of Q s . But the results of the fit do not change largely by changing Ao- The 
saturation exponent A is computable in QCD (known up to the renormalization-group- 
improved NLO BFKL) [TTJ [121 CHJ , but we treat A as a free parameter since the results 
are sensitive to its precise value. We work with three quarks of equal mass m q and use 
the same photon wavefunctions *$?t,l as in Refs. [3 [Hj- Thus, the only free parameters 
of the CGC fit are R, xq and A, which are the same as in the GBW model (ao = 2nR 2 ). 
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The CGC fit has been performed for the F 2 data at ZEUS [3] with x < 10~ 2 and Q 2 
between 0.045 and 45 GeV 2 (156 data points). In Fig. 1, the results of the fit are plotted 
against the data for Afo = 0.7 and m q = 140 MeV. The three parameters are determined 
to be R = 0.641 fm, x = 0.267 x 10~ 4 and A = 0.253 with rVd.o.f. = 0.81. We also 
show (with dashed line) the prediction of the BFKL calculation without saturation, as 
obtained by extrapolating the formula in the first line of Eq. (J2J) to arbitrarily large 
rQ s . This pure BFKL fit shows a too strong increase with 1/x at small Q 2 . On the 
other hand, the complete fit, including saturation, works remarkably well even at the 
lowest values of Q 2 that we have included. Note also that the value A = 0.25 determined 
from the fit is in good agreement with the theoretical result [13]. We have done the 
fit separately with the pure scaling part, and found that the fit becomes worse. This 
means that the "diffusion" term which violates the geometric scaling is crucial to fit the 
HERA data. This is not surprising because this term effectively changes the anomalous 
dimension from its BFKL value 7 S = 0.63 (for relatively large dipole sizes ^ 1/Q S ) 
to the DGLAP value 7 = 1 (for small dipole sizes), and thus partially simulates the 
DGLAP dynamics. However, as is evident from the figure, there is a deviation between 
the CGC fit and the data at high Q 2 and not so small x. This is again not surprising 
simply because this regime is outside the range of validity of the CGC fit, which does 
not include the DGLAP physics (in its right form) nor valence quark dynamics. 
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